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We introduce a new approach for studying unified supergravity models. In this approach 
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all the parameters of the grand unified theory (GUT) are fixed by imposing the corresponding 

Ph 

<D number of low energy observables. This determines the remaining particle spectrum whose de- 
pendence on the low energy observables can now be investigated. We also include some SUSY 
threshold corrections that have previously been neglected. In particular the SUSY thresh- 
old corrections to the fermion masses can have a significant impact on the Yukawa coupling 
unification. 



1. Introduction 



It has been shown recently that in the minimal supersymmetric model (MSSM) [1] the 
SU(3) C x SU(2) L xU(l) Y gauge couplings unify at a scale M GU T = C(10 16 GeV) [2]. Addi- 
tionally, the unification of r and bottom Yukawa couplings at Mgut can be achieved within 
the MSSM if the top- Yukawa coupling is close to the Landau-pole [3;4]. 

The evolution of the coupling constants from the electroweak scale, given by ra z , to the 
GUT scale, Mgut> including SUSY threshold corrections is non-trivial without any a priori 
knowledge of the SUSY particle spectrum. This particle spectrum results from mass parameters 
that are subject to a renormalization group (RG) evolution from Mgut to ra z and are in general 
scattered over several orders of magnitude. Several different approaches have been proposed, 
all of which make some simplified assumptions about the SUSY threshold effects. Also, it is in 
general not possible to impose all the experimental observables. For example, in the so-called 
bottom-up approach [5] which is suited to investigate a large number of models a s is an output 
rather than an experimental input. Any experimental information on a s is thus lost. 

In general, the goal of all these approaches was to obtain low energy predictions by impos- 
ing GUT constraints. In this paper we introduce a complementary approach that enables us to 
constrain the GUT parameter space by imposing all present and future experimental results. 
Furthermore, we present a more complete treatment of the SUSY threshold corrections that 
does not require any assumptions about the SUSY particle spectrum. Our approach is charac- 
terized by fixing all the low energy observables including the strong coupling constant, a s and 
the bottom mass, ra&, precisely. By varying these parameters over the range that is experimen- 
tally viable we still obtain the full range of viable SUSY parameters. This way we can explore 
the possibilities and limitations of probing the GUT parameter space by high-precision mea- 
surements. In turn, we can impose GUT constraints in order to study SUSY phenomenology 
in terms of only a few input parameters. 

Our strategy is as follows. We start out with a general N = 1 unified supergravity model. 
In the minimal SU(5) version there are nine input parameters at the GUT scale: the universal 
gauge coupling, a g , the up and down Yukawa couplings, ajj and ap, the Higgs mass parameter 
of the superpotential, /i, the GUT scale, Mgut> the universal mass parameter of the spin 
and spin 1/2 particles, tuq and m^ 2 , and the soft SUSY breaking parameters multiplying 
the trilinear and the quadratic part of the superpotential, A and B. (Actually, the Yukawa 
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couplings are 3x3 matrices, but we can neglect the couplings of the first two generation, which 
are small and whose origin is unclear. We will also not impose any constraints on the SUSY 
parameters coming from the non-observation of proton-decay since they are model-dependent.) 
In order to obtain the low energy parameters we will solve the RG equations (RGEs) twice. At 
the first run we will determine all the mass parameters and thus all the thresholds at which 
the RGEs have to be modified. At the second run we will evolve only the coupling constants 
under consideration of all the SUSY threshold corrections. This way we obtain all the low 
energy observables including those that are already determined by experiment. By imposing 
these experimental results we obtain a strongly constrained SUSY particle spectrum whose 
dependence on the few remaining parameters can now be studied. 

Our paper is organized as follows: in Section 2 we describe our approach with all the relevant 
corrections. In Section 3 we present our numerical results and in Section 4 we summarize our 
conclusions. The one-loop corrections to the fermion masses are given in the appendix. 

2. Derivation of the SUSY particle spectrum 

The derivation of the SUSY particle spectrum is straightforward. First we evolve all the 
parameters from Mgut to ra z using the SUSY (3 functions (two-loop for the couplings [6] and 
one-loop for the mass parameters [7]) where we impose universality on the mass parameters 
and unification of the gauge and Yukawa coupling constants. In order to assure that heavy 
particles with mass m decouple from the set of RGEs at a scale \fs < m we multiply the 
corresponding mass parameter in the (3 function by a step function, [0(x) = 0,1 for x < 
0,a: > 0, respectively]. Thus we set for the squark, slepton and Higgs mass parameters 

M? ^M?8(s-M?), (1) 



where i = Q n ,U n ,D n ,L n ,E n ,H m , [n = l,2,3;ra = 1,2), and we decouple the yl-parameters 
at scales below the average squark or slepton masses 

A Un - AuJ(2s -M\ -Ml), 

A Dn ^A D J(2s-Ml-Ml), (2) 
A En -+ A E J(2s -M) -M 2 ~). 
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For the Higgsino and gaugino mass parameters we have 



M -> M6(s - M 2 ) 



(3) 



where M = fi, M^, M g . With these mass parameters known, we can now compute the full 
SUSY mass spectrum. 

However, the coupling constants we have obtained sofar do not include any threshold correc- 
tions. These can be included by evolving the coupling constants a second time while changing 
the RGEs at every threshold. We start by running the gauge and Yukawa coupling constants, 
cti [i = 1, 2, 3, r, 6, t), from Mgut to the scale at which all strongly interacting SUSY particles 
decouple, Mqcd> using two-loop SUSY (3 functions. We define Mq CD = M^M^ because the 
decoupling of the top-squark has the strongest effect on the evolution of the Yukawa couplings 
and the quartic Higgs couplings which below Mqcd are allowed to evolve differently than the 
gauge coupling. The effects of the different squark and gluino thresholds on the gauge cou- 
plings, as well as the conversion from DR to MS scheme [8] can be easily incorporated at the 
one-loop level by writing 

(«r 1 )_ = («r 1 ) + -A l , (4) 

where the subscript + ( — ) denotes the value above (below) threshold, and we have 
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In the MSSM the Higgs potential is given in general by 



V = mi!*}*! + ml 2 §\§ 2 ~ [mi 2 $I$2 + h.c. 
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+ ±Ai($l$i) 2 + |A 2 ($J$ 2 ) 2 + A 3 ($ t 1 $i)($5$ 2 ) + A 4 ($ t 1 $ 2 )($J$ 1 ) 

+ {±A 5 ($ t 1 $ 2 ) 2 + [A 6 ($ t 1 $ 1 ) + A 7 ($5$ 2 )]$ t 1 $ 2 +h.c.} , (6) 

where the mass parameters are given at Mqut by = ra 2 ^ = mg + fi 2 and ra 2 2 = —fiB and 
have been evolved down to Mqcd- The coupling constants are given at Mqcd by 

Ai = \{g\ + g\) + AAi , 
A 2 = \{gl + fl'i) + AA 2 , 
As = \{g\ - g\) + AA 3 , 
A 4 = -iflr| + AA 4 . 

Here, the AA^ [i = 1, .., 7) are finite one-loop threshold corrections presented in ref. 9. At scales 
below Mqcd the quartic couplings of the Higgs potential will evolve according to the (3 function 
of a general two-Higgs-Doublet model (plus sleptons, charginos and neutralinos) [9;10]. With 
these (3 functions the couplings evolve down to the scale M weak = M^ Mg 2 where all the electro- 
weakly interacting superpartners decouple (M^ and Mg 2 denote the selectron mass eigenvalues). 
We also assume the m^o = 0(M wea ] !i ) so that only the standard model (SM) particle content 
remains at scales below Af^. Thus, for the evolution from M wea k down to ra z we use the SM 
(3 functions. The SM coupling constants are obtained by imposing the boundary conditions 

\ SM = cos 4 /3Ai + sin 4 (3\ 2 + 2 cos 2 (3 sin 2 (3\ z 
+ 4 cos 3 (3 sin (3\q + 4 cos (3 sin 3 f3\f , 
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where A3 = A3 + A4 + A5. In order to decouple the top Yukawa coupling from the set of RGEs 
for s < mt we write 

at — > at6(s — m? t ) . (9) 
Again we use eq. (4) to account for different thresholds, with 

40tt I \M 2 V J \M 2 J 

^ v weak ' \ weak / 



A5 = AA5 , 
A6 = AA6 , 
A7 = AA7 , 



(7) 
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(10) 



where at this point we approximate the CP-odd Higgs mass by ra 2 ^ = fn\\ + m 22 an( ^ the 
Higgsino mass by = fi. It is understood that all the masses lighter than ra z should be 
replaced by ra z . At the scale ra z we impose the minimum conditions for the Higgs potential of 
eq. (6) 



\v 2 cos z j3 ( Ai + 3A 6 tan j3 + A 3 tan z j3 + A 7 tan d j3 ) - tan /Sm^ = , 



sin 2 j3 ( A 2 + 3A 7 cot /3 + A 3 cot 2 /? + A 6 cot 3 /? + ra 2 ^ - cot (3m\ 2 = 



(11) 



where the low energy effective couplings are given by 



Ai = Ai(M weak ) + cos 4 /3AA 
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(12) 



with AA 5M = A 5M (m z ) - A 5M (M weak ). This way we have achieved that only the SM Higgs self 
coupling evolves at scales below m^o . By using running parameters in the Higgs sector we have 
included all leading log terms summed to all orders in perturbation theory. This formalism 
yields a better approximation for large SUSY masses [say 0(1 TeV)] than e.g. the one-loop 
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effective potential. The top quark threshold corrections to a s are included via eq. (4) with 



A3 = — In I — i 
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(13) 



Instead of solving eq. (11) for tan (3 and v we keep tan (3 fixed and solve for ra 2 2 and v. This is 
possible because until now we have not used the value of B which does not enter the (3 function 
of any other parameter at one-loop. It can thus be chosen independently of all the other 
parameters. We now obtain the experimental observables 

tan W = - — , 

"em = sin 2 W OL2 , 

a s = «3 , 
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where the superscripts t and z indicate that these are the running masses evaluated at mt and 
ra z , respectively. Included are also all the tan/3 enhanced one-loop SUSY threshold corrections, 
Amj and Am r , listed in the appendix [11; 12]. The top Yukawa coupling is evaluated at mt 
and all other couplings at ra z . The physical (pole) masses are [13] 
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We now have to compare out theoretical predictions with experimental data. Our set of exper- 
imental quantities is [14] 



sin 2 9 W = 0.2324 - 1.96 x 10 -3 



138 GeV 
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ra z = 91.187 GeV 



m T = 1.7841 GeV. 



(16) 



The experimental errors of these quantities are so small that they can be ignored. On the 
other hand, the QCD gauge coupling carries a significant error. We typically choose a s = 0.13, 
which is somewhat larger than the result from deep inelastic scattering, a s = 0.112 [15], or even 
from LEP experiments, a s = 0.122 [16]. This is necessary in order to obtain a SUSY particle 
spectrum at or below 1 TeV. This discrepancy which is still within the experimental bounds 
might also be an indication for non-negligible GUT threshold corrections [17]. The bottom 
mass has been calculated to be ra& = (4.72 ± 0.05) GeV [18]. We are more conservative and 
accept mj = 5.0 GeV as the upper limit still compatible with the experimental data. In order 
to avoid that the discrepancy in a s propagates into our prediction of ra& (which for a s = 0.13 
and rat/ sin /3 < 200 GeV is unacceptably large) we fix in our plots the running bottom mass 
at ra z , ra| = 3.3 GeV, which corresponds to a physical mass of ra& = 5.0 for a s = 0.115. 
Furthermore, we typically choose tan /3 = 35 and rat = 185 GeV (this corresponds to a running 
mass of about 172 GeV which is significantly below the IR fixed point [19]). 

Now we proceed to fix the GUT input parameters, denoted generically by 7^, by imposing 
the set of experimental observables described above and denoted by Oj. Sofar we have described 
how to obtain these experimental observables as functions of the GUT inputs, Oj = Fj(Ii). 
However, in practice we need to find Ii as a function of Oj. Thus, we have to invert Fj(Ii) 
numerically. Clearly, this requires the number of inputs to be equal to the number of outputs 
and can be achieved by defining Ii as the limit of an infinite series 



where we choose a set of initial GUT parameters, /?, by default. All the other elements of the 
series are obtained by iteration 



h = Fr\0 3 ) = lim I\ 
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(n = 0,l,... 



) where we have computed the Jacobian matrices 
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(19) 
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A / m [at M GUT ] 

Fig. 1 The low energy SUSY mass parameters as a function of the yl-parameter at Mgut- 
The low energy value of At is negative unless its initial value at Mgut 1S very large (> 1.2rao 
in this plot) 

and their inverses, numerically. Operationally, we truncate the series when the relative 

error, [Fj(I?)/Oj — l] 2 , is below the desired accuracy. Whether the series defined in eq. (18) 
converges or not depends strongly on the choice of the initial GUT parameters, 7°. Thus, the 
procedure might have to be repeated several times before a set of GUT parameters can be 
found that yields the desired low energy observables. 



3. Numerical Results 

In the last section we have seen how to obtain the GUT parameters and thus also the 
full SUSY particle spectrum by imposing low energy observables. We will now present the 
numerical results of our approach. 

3.1. Without GUT threshold corrections 

Let us assume for now, that the heavy GUT particle spectrum is roughly degenerated and 
hence that the GUT threshold corrections are negligible. We present our results as functions of 
mj and a s due to their large experimental errors and of vat and tan /?, which are still unknown. 

A few clarifications are required before presenting the numerical results. In the absence 
of large SUSY threshold corrections the value of ra& is predicted once we fix m T and vat- Let 
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m t [GeV] 

Fig. 2 The SUSY mass parameters as a function of the top mass. At is positive (negative) for 
m t > (<)193 GeV 

us now look at the one-loop contributions to ra& [eq. (14)] given in Appendix A. In the limit 
^> m w we obtain for the dominant gluino and chargino contributions 
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where m = max{M^ , , M g }. Thus, in the large tan (3 regime (which is under study here) 
the corrections are strongly enhanced and as a result the bottom mass prediction depends 
on various other SUSY parameters. Therefore, we can treat ra r , mi and mt as independent 
parameters. 

It has been shown, that the prediction of mi from unification without threshold corrections 
tends to be too large unless the top mass is close to its IR fixed point [3;4]. Thus, we will only 
obtain acceptable values of mi for fiM g < [in our plot we will choose M g > 0] . By looking at 
the RGEs we see that M g tends to drive At to negative values. As a result, there is a partial 
cancellation of the terms in eq. (20) over a large portion of the parameter space. Nonetheless, 
the corrections can be significant. 

In fig. 1 we present the SUSY parameters as a function of the yl-parameter at Mgut- We 
see that the SUSY particle spectrum changes only slowly as long as \A\ ^ mo. The reason it 
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Fig. 3 The SUSY mass parameters as a function of the running bottom mass, m^. The 
parameter At is positive (negative) for m^ > (<)3.45 GeV. 

that the evolution of the A parameters is dominated by M g and thus the low energy value is 
rather independent of the input at Mgut which we will choose to be A = mo. In fig. 2-5 we 
present various low energy SUSY parameters as a function of rat, ra&, tan (3 and a s , respectively. 
Shown is the region in parameter space where we can find a set of GUT parameters that yields 
the desired values of the experimental observables. We display the dependence of Mgut> which 
is important for proton decay, the gaugino and Higgsino mass parameters, M2 and /i, that 
determine the chargino and neutralino sector, the mass parameter of the left handed top and 
bottom squark, Mq , which characterizes the scale of the squark masses and determines the 
radiative corrections [20] to the lightest Higgs boson mass, ra/jo, and the mass of the CP-odd 
scalar 



m 2 A o 



2m\ 2 



\v 2 (2A 5 + A 6 cot j3 + X 7 tan j3) 



(21) 



sin 2/3 

We also present the 5-parameter, which turns out to be an important measure for the amount 
of fine-tuning required to yield large values of tan/3, and At = Ajj 3 for completeness. If one 
neglects the SUSY threshold corrections [eq. (20)] only a very narrow range of mt close to the 
IR fixed point is compatible with r-bottom Yukawa unification [3;4]. However, in fig. 2 we 
see that by including these corrections we can find a solution for a much larger range of top 
masses that is clearly below the IR-fixed point mt ~ 200 GeV [19]. It has been pointed out 
in ref. 11 that the corrections in eq. (20) are constrained by requiring the absence of severe 
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Fig. 4 The SUSY mass parameters as a function of tan (3 

fine-tuning. The reason is that the minimum conditions of the Higgs potential require that 
Bfi/m,Q = 0(1/ tan /3) <C 1. Such small values for Bfi are instable under radiative corrections 
and, therefore, unnatural unless they are protected by an approximate symmetry. Such a 
symmetry would also require At and m^ 2 to be suppressed which in turn would imply a 
cancellation of the tan (3 enhancement of eq. (20). However, our main priority is to look for 
solutions for a particular set of experimental inputs. Thus, we consider the amount of fine- 
tuning required for a particular solution as an output which is roughly characterized by the 
ratio r = maxjyl 2 , M$} / B 2 . The absence of fine-tuning is an important argument in order 
to find the theoretically favored range in parameter space. For example, the requirement that 
r ^ 100 implies that mt 180 GeV. If we ignore the problem of fine-tuning we find a lower 
limit of mt ^ 150 GeV by requiring that m^o > 0. It is interesting to note that the SUSY 
particle mass spectrum becomes heavier with mt as a result of the mt dependence of sin 2 # w 
[eq. (16)] and of a s [eq. (13)]. 

In fig. 3 we present the bottom mass range for which a solution can be found. It illustrates 
the significance of the radiative corrections to the r and b quark mass [Appendix A]. It is 
interesting that the ra& dependence shown in fig. 3 exhibits very similar qualitative features to 
the mt dependence in fig. 2. 

In fig. 4 we show the tan (3 dependence of the SUSY particle spectrum which is very similar 
to the mt and ra& dependence in fig. 2 and 3. The main difference is that the condition m^o > 
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Fig. 5 The SUSY mass parameters as a function of the strong coupling constant, cx. rms 

puts an upper limit on tan (3 ^ ht/hi, [21] as a result of radiative symmetry breaking [22] . 

In fig. 5 we display the range of the strong coupling constants for which gauge and r-bottom 
Yukawa unification can be achieved. We find that Mgut grows with a s while the SUSY mass 
spectrum decreases with increasing a s [2]. However, the actual SUSY particle masses turn out 
to be much larger than the effective SUSY scale [4]. As a result, the predicted values of a s from 
unification are significantly larger than the experimental value as long as the SUSY particle 
spectrum is not much heavier than 1 TeV. This might be a first indication that the GUT 
threshold corrections [17] may be significant. 

3.2. With GUT threshold corrections 

We now consider the possibility of non-negligible GUT threshold corrections to the gauge 
couplings. We parameterize these corrections by A such that 

CKj 1 = a 2 1 = a 3 1 — A = <Xg 1 . (22) 

Here, Mgut is defined as the point at which a\ and «2 intersect. It does in general not 
correspond to the masses of any heavy particles which are now assumed to be non-degenerate. 
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Fig. 6 The SUSY mass parameters as a function of the top mass. At is positive (negative) 
for mt > ( < ) 178 GeV. Included are finite GUT threshold corrections to the gauge couplings 
parameterized by Air A. Instead, we fix the lightest top squark mass, = 1 TeV 



In general, we can write 

4,A = Y / H^( 1 ^-) , (23) 

where the sum is over all heavy GUT particles, </>, and the coefficients = 0(1). For simplicity 
we set a.i = a T at Mgut- I n fig- 6 we present various low energy SUSY parameters as a function 
of mt. The difference to fig. 2 is that we now fix the lightest top squark mass, = 1 TeV, 
and instead we allow non-zero values of A. This allows us to fix the strong coupling constant at 
a lower value a s = 0.115 and to fix the physical (pole) bottom mass ra& = 5.0 GeV. The solid 
curve shows the predicted GUT threshold effect defined in eq. (23). We see that the full range 
of mt from the experimental lower bound to the IR fixed point is compatible with r-bottom 
Yukawa unification within the frame-work of the minimal unified SU(5) SUGRA model. 

3.3. SUSY LOOP-EFFECTS 

One of the main advantages of our top-down approach is that it allows us to greatly reduce 
the number of free SUSY parameters. Consider, e.g. the mass of the lightest Higgs boson, m^o. 
In the MSSM there exists a well defined tree-level upper limit m^o < m z which implies that 
the Higgs sector will be a good testing ground for the MSSM even if all SUSY partners are 
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Fig. 7 The two-loop radiatively corrected lightest Higgs mass (a) and some relevant SUSY 
particle tree-level masses (b) as a function of mt. We fix = 1 TeV 

heavy. Through radiative corrections m^o depends on all SUSY parameters [20] even though 
not all of them are significant. By imposing GUT relations, universality and the experimental 
constraints from eq. (16) the number of free parameters can be reduced to five. Furthermore, 
without great loss of generality we can set A = mo, since the low energy value of At is primarily 
determined by rrig. In fig. 7 (a) we present the two-loop radiatively corrected prediction of 
ra^o [20;23] as a function of mt. Here we set tan/3 = 35, Mj = 1 TeV and we have imposed the 
constraint fi = — mo/2. This constraint is chosen for convenience since it allows for radiative 
symmetry breaking over a large range of mt. In fig. 7 (b) we present some physical SUSY 
particle masses (at tree-level). We see that the constraints Mf x ,m^,m^± ^ m z /2 imply 
120 GeV ^ mt ^ 170 GeV. Of course, this range of mt for which a solution can be found 
depends on our initial values of A and fi. However, the value of m^o as a function will not be 
affected significantly. Here the QCD coupling and the bottom mass are treaded as outputs and 
lie in the range 0.122 < a s < 0.133 and 5.4 GeV < ra& < 6.6 GeV where the lower (upper) 
limits correspond to mt = 120 GeV (170 GeV). 
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4. Conclusions 



We have introduced a top-down approach for studying the GUT parameter space in SUSY- 
GUT models. With this approach we can analyze explicitly the dependence of the SUSY particle 
spectrum on various experimental observables. We find that the SUSY particle spectrum 
becomes heavier with increasing mt and decreasing a s . Neglecting GUT threshold corrections 
and assuming all SUSY particles are lighter than 1 TeV we find that a s ^ 0.13. In addition, we 
have included in our analysis the potentially large SUSY threshold corrections to the quark and 
lepton masses previously neglected. We focus our attention on the large tan (3 limit where we 
have shown that these corrections can in some cases lower the value of mt for which r-bottom 
Yukawa unification is achieved below the present experimental lower limit. 

Furthermore, our approach enables us to fix the entire SUSY particle spectrum in terms 
of a few experimental inputs. This allows the computation of virtual SUSY effects in terms of 
only a few parameters. 
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one-loop SUSY threshold corrections to ra& and m T 

In this appendix we present results for the one-loop radiatively generated down-type fermion 
masses. The calculation was done in n dimensions using dimensional regularization. We 
consider the limit of large tan/3 while keeping the Yukawa couplings h u oc m u / sin /3 and 
hd oc m^/cos/3 constant (here the index d stands for all down type fermions, in particular b 
and r). These types of corrections are finite and scheme-independent. The result for the up 
sfermion-chargino loops is 



n,i 

where n,z = 1,2. The vertices in the basis of electroweak eigenstates are given by [24] 

/ gm d Ui 2 \ 



APPENDIX 





V. 



mt 



L± 




m 
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= ( gmuV t2 J • (A.2) 
^ \/2m w sin (3 ' m 

The chargino mass eigenvalues, M-±, and rotation matrices, U and V, are denned via 



M m w V2 N 
-ft 



diag(M- i± ,M.±) = f/XF" 1 , X = ( W ) , (A3) 



From here we obtain the vertices for the mass eigenstates by rotating 



where P = L,R and the mixing angle and the unitary matrix are defined by 

2( Aj — cot f3)m v 



sin 26*5 = 5 2 ' 

to£ — to£ 

(cos sin 9 \ 
■ 
— sin 6 cos / 



(A5) 



Within the framework of a large tan (3 approximation we set = for consistency and 
finiteness. In this case, the conventionally defined scalar two-point function becomes 

„ , 9 o> . to? In to? — to? In to? . , . 

5 (0,to 2 ,to 2 ) = A-^ \ \ 2 - + 1 , (A6) 

where A = 2/(4 — n) — 73 + ln(47r) and 73 is the Euler constant. The result for the down 
sfermion-neutralino loops is 

i,n 

where i = 1,2,3,4 and n = 1,2. The vertices for the mass eigenstates are obtained by rotation 

• oa - 2(Ai -/^tan/3)TO d 

sin 20? = 5 5 , (A8) 

d to 2 - - to 2 - v ' 

dl 0,2 
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where we have defined 



V, 



V, 



LO 



x /2ee d Z' il -^-(l + 2e d s 2 w )Z' i2 

I c w 



i?o _ ^ [ m w cos /? 



72 



-2ee d Z\ x + 2-*-e d slZ< 2 



(A.9) 



Here the neutralino mass eigenvalues, -M^o, and rotation matrix, Z, are defined via 

diag^M^M^M^ZyZ" 1 , 
/ M' m z s w \ 

M -m z c w 

-^ 

\m z 5 w -m z c w -//, / 



(A.10) 



The result for the down squark-gluino loop is (these type of corrections are absent for the 
leptons) 



Am d = —Ms sin 26 ; 
d 3tt 9 d 



B (0,Ml,Ml)-B (0,Ml,M$) 



(All) 
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